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AN IMPROVED HELIUM-3 NEUTRON SPECTROMETER 



William A. Bair 

Lawrence Radiation Laboratory 
University of California 
Berkeley, California 



ABSTRACT 

3 

This report is a preliminary analysis of a He neutron spectrometer 

being developed at the Lawrence Radiation Laboratory, Berkeley. The 

spectrometer will be capable of presenting neutron spectra to an upper 

3 

limit of 14. 0 MeV. The spectrometer employs a He-Kr gas mixture in 

3 3 

a proportional counter, and the reactions He(n, p)T and He(n, d)D are 

3 

utilized to obtain the incident neutron energies. Recoils from He, Kr, 
and the stabilizing gas (CO 2 ) are electronically removed from the spectra 
by taking advantage of the shorter rise times of the recoil pulses. An 
analytic analysis of the rise time discrimination pulse rejection, the wall 
effects, and the spectrometer efficiency is presented. 

The computer program for the removal of the above corrections is 
lengthy, and although an outline of the methods of machine computation 
is presented, this concept is not carried to fruition. Instead a computer 
< program and an "inversion matrix" are developed (based on the analytic 
analysis presented) by means of a "random probability analysis." The 
"inversion matrix" should permit the correction of a raw neutron spectrum 
output from a multichannel analyzer to a final neutron spectrum by 



-Vll- 



correcting in a single operation for: (a) rise time discrimination, 

(b) wall effects, and (c) spectrometer efficiency. 

Appendices include a table of comparative physical and operational 
characteristics of all He proportional counters and ionization chambers 
constructed, and reported in the technical literature to date (Fall 1965). 
The complete electronic circuitry necessary to construct the rise time 
discriminator units, and the computer program "CATMAN" developed 
for the random probability analysis are also presented as separate 
appendice s. 
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I. INTRODUCTION 
A. Background and Purpose 

This report concerning an improved 3 He neutron spectrometer is 
a preliminary report concerning the theory, design, and construction 
of a ^He neutron spectrometer by the Health Physics Department of the 
Lawrence Radiation Laboratory. The actual operation and performance 
of the spectrometer will be the subject of a later report by the Health 
Physics Department and will be authored by Mr. Wai-Kit Quon. 

The detection of neutrons and especially the measurement of their 
energy have from the first attempts presented a more difficult problem 
than the other common nuclear particles. This difficulty is because the 
usual methods depend on effects resulting from charges on the particles. 
Most methods of neutron detection depend upon the scattering of particles 
after a collision with the neutron, and the measurement of the energy 
deposited by the scattered charged particles. The drawback to any 
method that is designed to measure neutron spectra which is based on 
scattering lies in the conversion of the definite energy of 'the incident 
neutron to a continuum of energies of the recoiling particle. To avoid 
this continuum of energies, a neutron induced nuclear transformation 
was sought in which the energy release would correspond unambiguously 
to the neutron energy. Investigations of several potential reactions, 
and the elimination of those that had sharp resonances, low lying energy 
states in the residual nucleus, and large positive- -or any magnitude 
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negative Q values, led to the selection of the following reaction: 

1 n Q + 3 He 2 — > 3 H 1 + 1 H 1 + 0. 764 MeV. 

In addition, the ^He(n, p)T cross section is fairly large and varies 
smoothly with energy, and the ability to discriminate against 7 rays 
enhances this selection for the task of neutron spectroscopy. Now, if 
the charged particles produced by the neutron interaction were to be 
completely stopped in the active volume of a proportional counter, the 
ionization caused by the loss of particle energy would be equal to the 
energies of the charged particles. This energy would then appear as 
free electrons and positive ions within the active volume. If an elec- 
trical field gradient is present, the electrons would migrate to the 
anode, and the charge deposited by the electrons would be proportional 
to the energy, of the incoming neutron. 

Several experimenters have examined the neutron interactions with 

the He nucleus, and have constructed neutron detectors, and to a 

1-13 

limited extent, neutron spectrometers utilizing proportional counters 

and ionization chambers. The work of the majority of the previous 

experimenters has, in general, concerned itself with the neutron energy 

spectra of only a few MeV, although Sayres reports success with neutron 

energies up to 8 . 1 MeV'*'® and the measurement of a cross section at 
5, 14 

17. 5 MeV using monoenergetic neutrons. 

In the majority of the previous work accomplished, two problems 
seem to inhibit the usefulness of a ^He neutron spectrometer. 
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3 

First, an ambiguity arises from the fact that He has a relatively 

3 3 

large cross section for elastic scattering. The elastic He(n, n) He 
cross section is equal to the nonelastic cross section at . 072 MeV, but 
increases rapidly with increasing energy. At 4. 2 MeV the elastic cross 
section is 4. 6 times as large as the inelastic cross section. The ratio 
decreases slightly with further increasing energy, but the elastic cross 
section is at least double the inelastic cross section from 0. 25 MeV to 
14 MeV. The probability of an undesired recoil pulse is thus at least 
twice the probability of a desired inelastic collision. Unfortunately, the 
recoil events, in addition to being more probable, also produce a 
continuum of energy pulses from 3/4 E n to zero depending upon the angle 
that the incoming neutrons strike a %e nucleus. This prevents utilizing 
any simple subtraction process in the determination of the initial neutron 
spectra. 

Second, the possibility of a nuclear event occurring in an orienta- 
tion such that the charged particles produced by the event intersect the 
walls of the containing vessel before expending their full energy was 
recognized. This phenomena gives rise to a count that is of less energy 
than a count that goes full track in the gas medium. He nonelastic 
collisions with high energy neutrons produce tracks that are of consider- 
able length in relation to the dimensions of the proportional counter, and 
the wall correction effects increase significantly with increasing neutron 
energies. 
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Solutions to the above two problems were attempted in various 
8, 9, 15 

manners. Brown constructed a "wall-less" counter composed 

of a peripheral ring of 1 6 small counters surrounding a central count- 
ing wire and an anticoincidence circuit to remove the wall effects. He 
then developed a computer program to unfold the He recoil events. 



16 



Brown 1 s work was extended by Wang who developed a computer program 



to unfold the "wall effects" in addition to the "recoil events. " Sayres 



10 



and Coppola developed a "risetime discriminator" circuitry that 



removes the recoil events {and a fraction of the desired events) by taking 



advantage of the shorter track length of a ^He recoil in comparison with 



the larger track lengths of the inelastic particles produced by neutron 



3 3, 4 

interaction with a He nucleus. Batchelor ' uses an approximate method 



of computation to remove the wall effects which appears to be useful at 
low energies, and short range of track in relationship to the counter 
dimensions. 



The purpose of this paper is to attempt to combine the removal of 



the %e(n,n)%e recoil effect developed by Sayres^ with the unfolding 



9 16 

of the wall effect developed by Brown and Wang. It is further the 



purpose of this paper to report attempts to extend the energy range of 
the He proportional counter to 14. 0 MeV (which appears to be about a 
practical upper limit for these devices) and to develop a means whereby 
the spectrometer can be utilized with some degree of 4 jt efficiency. The 



spectrometer will be used by the Health Physics Department, Lawrence 
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Radiation Laboratory, Berkeley for radiation surveys --hence the 

requirement for 4rt sensitivity. 

B. Description of the Spectrometer 

3 

The experimental apparatus for this He spectrometer consists 
of a proportional counter, electronic components, and a computer 
program. 

1 . Proportional counter 

The proportional counter is a cylindrical tube with an active length 
of 15 inches, a 4-inch outside diameter, and a wall thickness of 0. 035 
inch. The diameter of the anode wire is 0.003 inch. The proportional 
counter is designed with independently adjustable field tubes. Gas 

o 

filling is 10 atmospheres of krypton, 2 atmospheres of high purity J He 
(less than 10"^®% tritium), ^ and 1/2% CO 2 as a stabilizing gas. The 
tube was specially constructed by the Texas Nuclear Corporation for 
this spectrometer. A line diagram showing the tube and pertinent 
dimensions are shown in Fig. 1. 




Field tubes 

are independently 

adjustable 
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A discussion of the design parameters and the de sign, variable s 
is included in Appendix A. 

2. Electronic components 

The electronic components are designed to be able to distinguish 
between a pulse generated from a He recoil and a He inelastic event. 

To understand the operation of the electronics (a detailed description | 

of which appears in Appendix A), it is first necessary to understand 
what the electronics are to accomplish. Basically, the length of track 
of a He recoil (or any heavy particle recoil), is dependent upon the 
mass of the incident particle, the mass and charge of the recoil particle, 
the energy of the incident particle, the scattering angle, and the composi- 
tion of the absorbing medium. (An analysis of the collision dynamics 
involved are included in Appendix B. ) In general it can be stated that 
for equal recoil energies, a heavier particle, or a particle of the same 
mass but a higher atomic number, will traverse a shorter distance in 
the same stopping medium than a particle of less mass. In a propor- 
tional counter, a track r s radial component from the center wire determine 
the length of time that it takes to collect the charge, i.e. the pulse rise 
time. Since the track length of the recoils is always less than the track 
lengths of the disintegration particles (which are the desired events), it 
is possible to select a rise time that will effectively exclude all recoil 

'■'The arrangement also serves to eliminate recoils from the krypton gas 
and the CC> 2 . 
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events, yet retain a certain portion of the desired events which will 
have a longer rise time because of their longer track lengths. 
Unfortunately, the rise time criterion also eliminates many M good 
pulses” whose tracks are orientca with respect to the center wire in 
such a way that their rise times are within the discrimination criterion. 
The electronic circuitry developed in Appendix A is designed to 
accomplish this pulse rise time discrimination. 

3. Computer program 

The computer program is developed to accomplish the n unfolding ,f 

of the wall effects and to restore those desired pulses which were 

electronically eliminated by rise time discrimination. The wall effects 

1 6 

computer program follows the method outlined by Wang, but is not 
identical. Differences are primarily in the computation of the probability 
function P(E— > E^) which is the probability per unit energy that the 
reaction products of energy E will deposit energy'E^ inside the sensitive 
volume due to the wall effect. In addition, a new probability function- is 
analytically developed to replace mathematically the valid pulses which 
were eliminated due to rise time discrimination. An inversion matrix 
is proposed which will combine the two probability functions, and enable 
a complete reconstruction of the neutron spectra, above a minimum 
energy that is determined by our rise time discrimator setting. The 
neutron spectrum below this minimum energy cutoff level is unfortunately 
lost in the process of recoil elimination. 
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Unfortunately, the computer program analytically developed in 
this paper requires so much computer time that it was not considered 
warranted to carry it to completion. Instead, a random probability 
analysis was carried out, and the inversion matrix was developed on 
the basis of these random probability calculations. An outline of the 
random probability method is included in Section II-E and Appendix E. 
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II. DEVELOPMENT OF THE 3 He SPECTROMETER 

A. Helium-3 and Neutron Interactions 

A neutron can interact with a ^He nucleus in five possible ways. 

These reactions which are energy dependent are as follows: 

3 

1. He(n, p)T for neutrons of all energies 

3 3 

2. He(n,n) He for neutrons of all energies 

3 

3. He(n, d)D for neutrons of energies 4. 36 MeV 

3 

4. He(n, pn)D for neutrons of energies -r? 7. 32 MeV 

3 l 

5. He(n, p2n) H for neutrons of energies > 10. 3 MeV 

The cross sections for these reactions are reasonably well known 

1 8 

in the energy range 0. 001 eV to 14. 0 MeV with the exception of the 

o 

^Hefa, pn)D for which no data are available. It is assumed that this 

exceptional case has a small cross section, even though it is energetic 

cally possible. A comparative chart of the other reaction cross sections 

is included as Fig. 2. In addition to the He(n, pn)D reaction, the 
3 1 

He (n, p2n) H reaction will not be considered in the determination of' 
either wall effects or rise-time discrimination. The reasons for this 
omission are: 

a. At the maximum energy expected to be encountered in this proportional 
counter (i. e. 14. 0 MeV), the (n, 2n) cross section is only 15% of the total 
non-elastic cross section . 

3 

b. The neutron induced disintegration of the He nucleus involves four 
separate product particle s - -only two of which are charged. Four particle 



; /W 
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dynamics are complicated, and the energy imparted to the various 
particles may span a wide range of variables. Note, however, that the 
threshold for this reaction is 10. 3 MeV and the kinetic energy of the 
charged particles produced by even a 14.0 MeV incoming neutron cannot 
exceed 3. 7 MeV. Since some of the energy of the incoming neutron will 
undoubtedly be carried away as kinetic energy of the two product neutrons, 
it is highly probable that the sum of the charged-particle kinetic energies 
will be less than 3. 7 MeV. Pulses from a 14. 0 MeV neutron (n, 2n) 
reaction would thus appear on a multichannel analyzer at less than 3. 7 
MeV. Referring to Fig. 2, it is to be noted that the cross section of the 
^He(n, p)T reaction is fairly constant between about 0. 4 and 2. 5 MeV, and 
, the cross section at these energies is at least 14 times as great as the 
(n, 2n) reaction at 14. 0 MeV. Distortions of the neutron spectra would 
thus appear above 10. 3 MeV at essentially 0% at threshold, up to a 
maximum of 15% at 14. 0 MeV. This distortion would also manifest 
itself as a lower energy pulse, (^3. 7 MeV) but the allocation of reaction 
particle kinetic energies would result in an energy distribution that would 
not be unique, but rather span from 0 to 3. 7 MeV. For a continuous 
neutron spectra this "lower energy distortion*' is not considered serious 
due to the more probable (n, p) reaction at these lower energies which 
would essentially mask these degenerate pulses from higher energies. 

If we were to consider only single interactions within the active 

■ 

volume of our container (i. e. no multiple . scattering events), and to 



(borns) 
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Energy (MeV) 



MU. 37193 



3 

Fig. 2. Energy Dependent He Neutron Cross Sections 
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likewise consider no neutron disintegration interactions with the material 
in the container walls, we would find that (after recoils are eliminated 
by rise -time discrimination and wall effects mathematically resolved)' 
three reactions which would still tend to disturb our spectrum (irrespecti\ 
of all corrections). These reactions therefore limit the maximum 
efficiency of our spectrometer. The three interactions are: 

1. The disturbance of the energy spectrum above 10. 3 MeV and below 

3 1 

3. 7 MeV from the He(n, p2n) H reaction. 

2. The disturbances of the energy spectrum above 7. 32 MeV and below 

6. 68 MeV from the He(n,pn)D reaction (if this reaction does in fact occu: 

3. The disturbance of the energy spectrum at all energies due to the 

* 

neutron induced krypton disintegrations. 

The cumulative error introduced by neglecting these reactions is 
unknown. 

The remainder of this paper is to be divided into four parts. The 
plan of "unfolding” the raw spectra to develop a complete neutron spectrun 
is as follows : 

Development of a probability function to replace the pulses eliminated 
by ’’rise-time discrimination. ” (II-B) 

Development of a probability function to remove the wall and end 

t'fi, 

This problem is discussed in detail in Appendix A, although no attempt 
is made to correct the spectrum. 
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effects. (II-C) 

Combination of the two correction factors to reconstruct the neutron 
spectra. (II-D) 

Reasons for and development of a random probability inversion 

matrix in lieu of the above II-D. (II-E) 

Finally, it should be mentioned here that the work leading to the 

mathematical expressions for wall and end effects is primarily a restate- 

9 i / 

ment of the work of Dr. W. R. Brown and Mr. H. T. Wang , although 

some modifications have been made to the analysis of the later author. 

A number of their drawings, graphs, and formulas are repeated in this 

text for continuity and to facilitate explanations. They will be designated 

by a "dagger 11 (y) if. attributed to Dr. Brown, and by a "double dagger" 

(yf)if developed by Mr. Wang. 

B. Development of a Probability Function to Replace Pulses 
Eliminated by "Rise-time Discrimination 11 

•As a result of the arbitrary orientations of the disintegration particles 

3 3 

of the He{n, p)T and the He(n, d)D reactions, some tracks will be posi- 
tioned such that the time difference between the arrival of the first 
electron and the last electron to the central counting wire will be quite 
short. This situation will occur when the radial component of the track 
lengths is small, or when the event takes place near the central portion 
of the counter where the electrical field gradient is high. In both of these 
orientations the count produced by a "valid event" will be rejected by the 
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electronic .circuitry as though it were a recoil pulse. It was considered 
that if an analytical expression could be developed for determining the 
time difference between the arrival of the first and the last electron to 
the central counting wire, the "rise time” of the pulses could be deter- 
mined mathematically. A probability. function could then be developed 
based on the registered counts that would be representative of those good 
events that were rejected with the recoils. To develop this probability 
function, it is necessary however, to make four assumptions: 

(1.) The time required for the charged particles to deposit their 
energy in the gas medium is instantaneous. 

(2.) A nuclear event can occur with equal probability within any unit 
volume in the sensitive region of the proportional counter. 

(3.) The mobility of the electrons in the gas .mixture is^ directly’ propor- 
tional to the field gradient over the limited range of the electrical field 
gradients in which we are interested. 

2 

(4.) The electron mobility in the gas mixture (2 atmospheres He, 

10 atmospheres Kr, and 0. 5% CO 2 ) can be computed by assuming that the 



'■'As an approximation, it can be shown that by relating a specific particle 
energy Eq = -r rmq ; solving for v^; and then relating an average velocity 
to track length (j?-(E)) and solving for time t^ = 



that the time required for a particle to deposit its energy is about 3 nano- 
seconds for a proton (longest track length, highest velocity) of 10. 0 MeV 
energy, and 0. 6 nanoseconds for a 10. 0 MeV ^He recoil particle (shortest 
track length, smallest velocity). Rise times of interest to our counter 
circuitry are in excess of 3 microseconds. 
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electrons behave as ions, and that the mixture mobility can be computed 
as can the ’’mixture mobilities” of ions. 

Assumptions (3) and (4) will be di ■ tssed in greater detail in the 
following sections. 

1 . Electron migration time from any point r in coaxial cylindrical 

geometry . In a cylindrical proportional counter, the electrical field 

gradient can be determined at any point r from the center of the counter 

by the relationship: y 

E(r) = 2 

„ r c 
rJin — 

r a 

where E(r) is the field gradient at any radial distance r from the 
anode. 



( 1 ) 



V = applied voltage 

r = radius of the cathode 

c 

r = radius of the anode 

a 

(This formula is developed in Appendix D. ) 

For our counter dimensions we can plot the field gradient at any, 
radius from the center wire by using the appropriate applied voltage and 
counter dimensions. Instead of directly plotting the field gradient as 
a function of radius however, let us divide both sides of Eq. 1 by the 
total pressure of the counter filling gas in cm of mercury and plot this 
value instead. The reasons for this will become evident shortly. This 
plot can be seen' in Fig. 3 for representative counter voltages. 

An electron produced by ionization at a radial position r will drift 
slowly towards the anode under the influence of the electrical field gradient. 



(volts /cm • cm -1 Hg) 
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mu. 37194 



Fig. 3. Variation of Field Gradient/Filling Pressure (E/P) 
for Representative Applied 'Voltage 
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As can be seen in Fig. 3, the electrical field gradient at regions 
beyond about 1 cm from the cathode varies very slowly. After the 
electron passes this point in its inward radial drift, the field gradient 
sharply increases. If we assume a direct proportionality between the 
drift speed and the *E/P ratio, we can obtain an expression in cylindrical 
geometry for the time it requires for an electron to reach the cathode, i. e. 



where m, the electron mobility is the slope of the speed vs E/P curve. 



that within the limited range of E/P which is the most important in deter- 
mining the drift time of electrons for cylindrical proportional counters, 
that this assumption of linearity does not appear to be unreasonable. 
Substituting the expression for field strength in the case of cylindrical 
coaxial geometry, we obtain the electron drift time from any radial 



r 




( 2 ) 



19 



Measurements 



of drift speed for ^He and Kr with 0. 5% CO^ indicate 



position r. 




r dr 




(3) 



2 



t 




Since _iL is very small, the time. necessary to travel from the surface of 
the anode to the geometrical center is negligible. 
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The development of Eq. 2 and 3 follows the logic developed by 
Harling. ^ 9 

We still must determine however, some means of combining the 

mobilities of the electrons in a n gas mixture 11 since data is available 

only for mixtures of individual noble gases and moderating gases. A 

theoretical development of this "gas mixture" mobility problem was 

2 1 

presented by Staub, however the necessary gas constants are not 
available to compute the final result. We must therefore make use of 
assumption (d) and assume that as a first approximation, the electron 
mobilities behave as do ionic mobilities. 

While it is realized that no simple theory can account for the electron 
diffusion, attachment, recombination, and agitation velocity, as the 
electron moves through a changing field gradient, it is considered that 
this approximation will cause at most an error of a constant factor. The 
electron mobilities are thus assumed to follow the same physical laws 
as do ions --although it is realized that the situation is considerably more 
complex for the fast moving electrons. 

2 1 

Staub 1 s method of combining ion mobilities is as follows: 

i _ - p i i +£2 . _i_ 

Pl2 P /*1 P pz 

where 2 = the ion mobility in the gas mixture 

^12“ the ion mobility in the component gases 

]? 1 , 2 = partial pressures of the component gases 
P = + P 2 = total gas pressure. 
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Substituting m (the electron mobility) for p we have: 



1 _ P He J_ P Kr 

m mix p tot m He p tot 



1 

m Kr 



( 4 ) 



Using the sources explained in the footnote, we find that for our tube 



m 



He 



m. 



Kr 



He 



Kr 



V _i 

= 0.3 cm/^-sec- (— - • cm Hg) 



V 



= 1.0 cm/p sec- ( ^^crrT 1 Hg) 

= 2. 01 atmosphere 

= 10.05 atmospheres 



Substitution of the above values into Eq. 4 yields m . = 0. 72 

mix 

, V - j 

crr^u-sec- ( cm 'em Hg). Utilizing this value in Eq. 3 and the appropriate 
counter dimensions we can plot migration time for an electron from any 
radial point in the counter. This has been accomplished in Fig. 4. for 
representative counter operating voltages. 

2 # Rise time. Rise time of the pulse has been defined as the time 
difference between the arrival of the first electron at the central counting 
wire and the arrival of the last electron. By use of assumptions (c) and (d) 



20 

'‘'The rrlobility factor for Kr was taken from the data of English and Hanna. 

The .mobility factor for the helium is taken from the data of Friedes and 

Chrien, ^ which is plotted in reference 19 (Harling). The reason for this 

use of Friedes and Chrien* s data rather than Harling 1 s is that Harling, 

(1) made his measurements using ^He with only small amounts of %e, and 

(2) he used substantial percentages of CC> 2 . Although drift velocity of 
electrons in ^He appears to be rather insensitive to moderating gases, the 
CC >2 content of Friedes and Chrien was considerably less than that utilized 
by Harling. 



Electron migration time (//.sec) 




r (cm) 



mu. 37196 



Fig. 4. Variation of Electron Migration Time from Various 
Radial Positions for Representative Voltages 



we have been able to determine migration time, and have formulated 



an expression for the time for an electron to transit any radial distance. 
Thus the rise time (t ) can be expressed as 



Note that this is identical to the time that it takes an electron at 
position r r (furthest distance) to migrate to r n (nearest distance) if we 



previously noted that the field gradient charges very slowly beyond about 
1 cm from the anode wire (see Fig. 3). Therefore, if we select our rise 
time discriminator in such a manner that pulses which spend their full 
energy within the active volume enclosed by a 1 cm radius. are rejected. 
We can see that the assumption regarding the linear dependence of 
mobility and E/P is not too serious a concept since the pulses which 
depart from this assumed linear dependence(at high E/P values) are 
automatically eliminated. As it turns out, the "dead volume" of our' 
counter will extend to about 2. 75 cm, as will be shown later. Pulses 
which originate in the "dead volume" and extend outwards will be treated 
in detail later, as will pulses which do the reverse. (Section II-B. 7) 

3.. Electron production. From the above, we now have a means of 
calculating rise times for electrons that are located at various radial 
distances in the counter active volume. We do not, however, have any 




charge the limits of Eq. 3 from r to r„ and r to r„. We have 

fa 11 



indication of where these electrons will be located--in fact there will be 



? f/L 
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\ 



millions of electrons generated along the path length of a charged 

particle that is stopped in the gas medium. Let us now turn our attention 

to the formation of these electrons in a gas medium. We have previously 

discussed the various reactions that can occur with a He neutron inter « 

action, and have determined that only three of the five potential reactions 

are of interest in our counter. Neglecting the energy distribution within 

the three events for the moment, let us now determine the ranges of 

the various particles in the gas medium. 

2 2 

Friedlander and Kennedy present an empirical relationship for the 

determination of the ranges of charged particles in a gas mixture. This 

relationship, adapted for our specific usage, can be used to determine 

3 

the ranges of protons, tritons, deuterons, and He recoils as follows: 



Range formulas for the individual gases (mg/ cm 
Protons 



2 



) 



Rp in Kr (E) = 

R p in He (E) = 

^He Recoils 

R 3 in Kr (E) = 

R 3 He in He 



R p in air (E) [l. 89 - 0. 25 log 10 (E)] - 0. 36 
Rp in air (E) [o. 82 + 0. 043 log^Q (E)] 



(6 



| jRp ^ air [l. 89 - 0. 25 log 10 (|j]- 0. 36 

1 Rp in air (f)[o.82 + 0. 043 log 10 ('£)] 



The energy distribution is discussed in Appendix C and in Section II-B. 6. 
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Tritons 

R t in Kr (E) 
Rrp in He (E) 

Deutorons 

R d in Kr (E) 
R^ in He (E) 



3 |r p in air (f-j [l . 89 - 0. 25 log 10 
3 Rp in air ( P-j [o . P2 + 0. 043 logp-j 

2 |r p in air (^) ■ 89 “ °- 25 lo §10 

2 R p in air ( [o. 82 + 0. 043 logjg 





The particle ranges can be converted into centimeters by use of the 
relationship: 



JL (cm) = 



R tot (mg/cm 2 ) 

KZjeEi) p] 



14. 7, 



14. 7 



Kr 



(7) 



where is the particle range in centimeters in the gas mixture, 
P = partial pressure of the individual gases (psi) 

P = gas density at STP 

= total particle range in the gas mixture (see below) 

_J _ W He + W Kr 

R tot< E > R He( E > R Kr< E > 

W = weight fraction, W Pe + W^ r = 1 



W 



He' 



W 



Kr = 



P He PHe 

P He P He + P Kr P Kr 
P Kr PKr 



( 8 ) 



PHe p He + ?Kr P Kr 



1 - W He 
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where 

P = partial pressure of gas (psi) 

P = gas density at STP 
'°He = Q-tbbiy g/liter (f) 

P Xr = 3.455 g/liter (t) 

The above relationships enable us to conveniently compute the 
ranges of the charged particle in the gas mixtures, and are amenable 
to machine computation. The ranges computed as above are plotted in 
Fig. 5 for comparison. Note the substantial difference in the length of 
the track for various particles of the same energy- -particularly above 
about 1. 0 MeV. It is this difference in track length that will enable 
us to remove the recoil spectrum from our spectrometer. 

4, Rise time discriminator values . The rise time discriminator setting 
for the pulses can now be determined since we now have a means of 
determining the pulse rise times and the length of track of various 
particles in the gas medium. If we recognize that it is only the radial 
component of the track length that is important in our determination of 
rise time de scrimination, we can select the mo st unfavor able recoil 
situation that may be encountered. Setting the discriminator at the rise 
time appropriate to this most unfavorable case, we will remove from 
consideration all pulses that can be attributed to recoils (%e, C, O, and 
Xr). This worst case can be considered to be composed of two factors; 

(1) The maximum energy that can be imparted to a recoil nucleus, and 

(2) the longest rise time that this maximum energy recoil can cause, due 



Length of track (cm) 
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Energy of particle (MeV) 

mu. 37197 



Fig. 5. Ranges of Various Charged Particles in Gas 
Mixture of 10 Atmos. Kr, 2 Atmos. ^He, and 
0. 5% C0 2 
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to its geometrical position in the counter. 

Appendix B (Case 1) summarizes the maximum energy that can be 
imparted to the various recoil nuclei. It can be readily seen that the 
3 He nucleus can absorb up to a maximum of 3 / 4 the energy of the incidei 
neutron. The other possibilities for energy absorption are of consider- 
ably smaller magnitude. The ^He recoil is thus selected as the recoil 

of concern. The ranges of the other recoil particles in the gas mixture 

3 

is likewise less for equal energies when compared with He. The worst 
geometrical case that can occur is if the track length of the recoil ^He 
nucleus is oriented in a directly radial position with one end of the 
track at an infinitesimal distance from the wall of the cathode. This 



positioning is illustrated in Fig. 6 . 

Rise time discriminator setting can then be determined by: 

( 1 ) calculating the maximum recoil energy that can be imparted 



( max ) = 0.75 (max) 



( 2 ) determining the path length in the gas medium, ij(E 3 r (max))from 



He 



the equations presented in Part II-By,2. 

(3) determine maximum rise time by substituting into Eq. 5 the value 
b, the radius of the cathode, for r£, and [b -^(E^ (max ))] for r n . 



This rise time discriminator setting (t ) will then exclude all recoils 



from the resultant spectra. 
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Fig. 6. He recoil track position to give maximum pulse 
rise time. 

5 % Determination of radial positions of the nearest (r n ) and furthest (rf) 

• 3 3 

electrons from the anode for Ke(n, p)T and He(n, d)D reactions. 

Consider the fact that to define the positioning of any event in space, it is 
necessary to define six degrees of freedom. Three of these degrees of 
freedom are necessary to define the positioning of the center of gravity, 
and three are required to explain the orientation of the event with respect 
to the center of gravity. For the purposes of discussion, let us assume 
that we can define the "center of gravity" as the position where the neutron 
strikes a ^He atom. Since the reaction occurs within the boundaries of 
the counting tube, we can define the positioning of the event (center of 
gravity) in cylindrical coordinates of radius, length, and angular dependence. 
The operation of the counter, however, is independent of the positioning 
of the event in terms of length (except for end effects, which will be 
discussed later) and completely independent of any angular relationship, 
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wit h respect to the anode wire. This means that the positioning of the 
center of gravity can be defined in terms of a single parameter r. 

To define the orientation of the resultant ionization tracks produced 
by °Ke disintegration let us first turn our attention to the He(n, p)T 
reaction. For a neutron of given energy E n , and given resultant proton 
energy Ep, we can determine the energy E- of the triton particle and the I 
laboratory angle C£j between the proton vector and the triton vector by 
use of the relationships developed in Appendix B. 



= E 



N 



Q 



n, p 



- E, 



and 



(2 E t + Q n p ) 
cosa, — 

^ 2 / 3 E t E p 



(1C 



Note that a is always obtuse in the laboratory system of coordinates. 

Let us now define a coordinate system. The coordinate system for 
our calculations is selected as follows: 

(1) The event occurs at a perpendicular radial distance rg from the 
center line of the anode. The end point of rg defines the zero point 
(0, 0, 0) of the coordinate axis. 

(2) The z axis is selected to 'include the vector r and extensions thereof. 

(3) The x axis extends parallel to the anode wire and passes through 



the point (0, 0, 0). 
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(4) The y axis passes through the point (0, 0, 0) defined as above, 
and is mutually perpendicular to the x and z axis. 



z 




Fig. 7. .Selection of coordinate system. 

With this selection of coordinate axis,' the event is angularly independent 
of the center anode wire. Calculations of the path orientations of the 
ionizing tracks, will be made in terms of polar coordinates using these 
axes as reference lines. We must now determine the distances of the 
nearest and farthest electrons from the anode wire. For the purpose 
of illustration let us visualize the Ep and Erp track lengths as solid lines. 
The orientation of these two vectors can have three degrees of freedom. 
Let us select the three rotational axes as shown in Fig. 8. 

Neglecting for the moment the E, p vector, let us turn our attention 
to the Ep vector. 

We are able to define a position which we call the tip of the Ep 
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Fig. 8. Rotational axes. 



vector. (We will refer to the Ep vector, although in reality we are 
speaking of a definite* length of path of a proton of set energy in the 
counter gas medium. ) We can completely describe this point r^ with 
respect to our coordinate axis by the polar angle 62, the asmuthal angle 
$2 and the length of the. vector Ep. We can thus determine ; r £ (Ep, 62* $2) 




Fig. 9 . r 2 (E p , Q 2 , 0 2 ) , 

We need an expression of r 2 in terms of rQ since in determining migra- 



tion times we are concerned with the radial distance from the central 
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Fig. 10. Determination of ^(Tq, E^, 9^, 0£) 
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counting wire to the point (Ep, 9^, 0 2 ), i. e. we need an expression 
of r 2 (r Q , Ep, 0^, 02)- To. determine this expres sion let us pass a. plane 
which contains the x axis and the vector Ep (a plane can be defined by 
two intersecting lines). This plane will intersect the z axis at (0, 0, 0) 
and will make an angle (3 with the z axis in the minus y direction. 
Referring to Fig. 10 we can derive an expression for ^(^q, Ep, 9^, 0^). 
3y the law of cosines: 

r 2 = r Q + ~ sin^ 9^ cos^(- 0 ^ )j " ^ r Q " sin^ 9~> cos^(- 0^)] ^ • cos (1 

but, cos (180 - ) = - cos (3 

Ep cos 02 

and from Fig. 10, - cos £= — — r 77/2 > 

Ep [JL “ sin^ ©2 cos^ (- 0 2 )j 

substituting this value for cos £3 into the above equation, we find that 



2 2 2 r 2 

= r + E [l - sin^ 0 O cos (- 0 9 )"| t 2 r A E— cos 9, 

2 0 P ^ ^ J OP L 



or 



(1 



r 2 ( r 0 » E p» Q Z’ ^2) = /^ Q + E p[ 1 ■ sin 2 9 2 cos Z (-0 2 )J + 2r Q E p cos 9 2 



Let us now consider the possibility that since we have not placed a 
restriction on where our event can occur, that a skew line may develop 
and that the closest point to the anode wire may in fact lie somewhere 
along the length of vector Ep, rather than at the end point. Since Ep 
is a definite length in the gas medium for a given energy of the proton, 
let us call this length J2(Ep) . An axiom of solid geometry is that the 
closest point of approach for skew lines is at a point where a line 
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connecting the skew lines is perpendicular to each. This means that 

a minimum distance can be determined. We have an expression for the 

length r^. Therefore, if we take a partial derivative of with respect 

to Ep (note that this is actually ^(Ep))and set it equal to zero, we can 

determine a di stance JLp along the X(Ep) vector whereby the skew lines 

are closest. We will call the closest distance between the vector Ep 

and the anode wire r ,. 

4 



o r- 



S4£to) 



“ ^r 2 -^(Ep) 2 [1 - sin 2 0 2 cos^ (- 02)] + 2i?(E p ) r Q cos 0 2 j 



fl/2 



' '2i(Ep) |l - sin 2 02 cos 2 (- 02 )j +2rg cos 0^ 



Set the above expression equal to zero and solve for Xp; note that the 
first bracketed term can never equal zero. For real values of Xp we 



o-otam 



= - 



r 0 cos 0 2 



P [l - sin^ ©2 cos^ (- 0 ^ )j 



( 12 ) 



where 0 < X p <i(Ep). 



The negative sign indicates that j£p cannot be positive unless ©2 is 
greater than 90°. This is in accordance with the physical situation 
envisioned and the selection of the coordinate axes. At ©2 = re, JLp must 
equal Tq and the Ep vector is perpendicular to the anode wire. Xp cannot 
exceed the physical length of Ep which has been assigned. 
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Substitution of this value of £ p into Eq. 11 for Ep determines 



that: 



r 4 = 



or 



2 r 0 cos2 0 2 • L 1 - sin2 e 2 cos2 (- ^2 )] - 2 r 0 cos2 0 2 

fr + • 



0 [jL- - sir? co (-02)/ 



[l - sin^ 0^ cos Z (- 0 2 )] 



r o 



1 - 



cos 9- 



1.1 - sin 2 02 cos 2 (- 0 2 ) j 



where — -^0 7 ^ " 

2 2 



and 0 ^ j2-(Ep) 



Let us now turn our attention to the other portion of the He (n, p)T 
event. Let us consider the Ep vector. 

In a manner similar to the determination of the end of the proton 
vector Ep,we can position the end of the triton vector Ep. We shall 




-X 



In the case of the triton however, we cannot arbitrarily describe the 
position of the end of the Ep vector with respect to the anode wire, since 



there exists a unique relationship for the positioning of the triton 
vector for a given neutron and proton energy. (Ejq and Ep respec- 
tively. ) This relationship is characterized by the angle OL which is 
the obtuse angle between the proton and triton. We must therefore 
consider this constraint upon the direction that the triton can travel. 
We can visualize this restriction by the following: 




Fig. 12. Rotation of Ej, Vector about Ep axis. 

If we were to rotate the Ep vector in Fig. 12 about its own axis 
while imagining that the E vector is rigidly attached to the Ep vector 
at an angle from the axis of the Ep vector, we can see that the tip 
of the Erp vector would trace a circle in a plane that is perpendicular 
to the Ep vector. The determination of the radial distance to the center 
wire r^(rg, ®1» ^0 mu take this degree of freedom into account. 

Let us now examine Fig. 13. The interior angie of the cone described 



-36- 



by the rotation of E.p is (180 - C^). Let us now pass a plane through the 
z axis and containing the Ep vector. Since this plane contains the z axi 
it will by definition be perpendicular to the plane containing the x and y 
axis. We will use this plane to relate 0^ to 9^, and 0^ to 0^. Let us 
consider the rotation of the E^ vector from a 11 zero position 11 where the 
Erp vector is at a maximum position away from the center wire (i. e. 
top position) and is contained in the (z, Ep) plane. The Ep vector is 
then rotated in a counter clockwise manner (viewed from the tip of the 
Ep vector). The angle of rotation we call 7 $ . By projecting the 
Erp vector back to the plane containing Ep and .z we can determine the 
incremental angle between 0^ and (180 - 02 ) this incremental angle is 



sm 



E sin a L cos y 
|E T (1 - sin 2 sin z 7 ) l/2 j 



. -1 

or sm 



sinaL cos 7 
(1 - sin 2 sin^ 



therefore 9^ is related to 9^ by 



i [ sin QIt cos y 

9i = (180 - 9 2 ) ± sin" 1 / ■■■■ ry 

1 2 'p - sin 2 Ci L sin 2 7 ) irZ 



or 



9 X = 180 - Q 2 + sin" 1 



sin <2 l cos 7 



,(1 - Sin 2 a L sin 2 7 ) 1/2 



(1 



where as determined by E^ T and Ep is a constant, 7 varies from 

0to2^. The sign of the incremental angleis positive for - JS < 7 <?L and 

2 2 

negative for 7 0 9 is arbitrarily selected. 

2 2 * 
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z 




Note: Vector E T shown exaggerated in size with 

respect to E p for clarity of angular relationships. 



mu. 37201 



Fig. 13. Determining Relationship Between 02 and. Op and 
02 and 0^ 
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If we likewise were to consider the perpendicular projection of the 
tip of the E rr. vector to the (Ep, z) plane we can readily determine that 
since this projection is perpendicular to the (Ep, z) plane, it must of 
necessity be parallel to x - y plane. We may therefore compute the 
incremental angle between 0^ and 0^, and determine a relationship 
between them. The incremental angle is 

f _ / i .? -2 v 1 / 2 •} 

i E'-n (1 - sin^C^T sm y) 

cos" 1 ± ±± — 

L j 

JL -j 

Therefore, the relationship between 0^ and 0^ can be expressed as 

(180 - 0^) = (180-0^) - cos"^{jl- sin^ CXj^ sin^ y)f ^ 
or 

0 j, = 02 i. cos * £(1 - sin^a ^ sin^y ^ J 

where the negative sign of the square root is taken for 0 \ y < and the 
positive square root is taken for it </ <2:t 

is a constant determined by fixed Ep and E^ 

02 is arbitrary 

Having now established a relationship between 0} and 02, 0} and 02 
through the constant angle c{-^and the rotational, angle y, let us now 
proceed to determine r ^ (r q, E,p, 0p 0^ ). 

Referring to Fig. 14, the same procedure can be followed as was 
used for the computation of r 2 (r Q , Ep, 0 2 , 0 2 ) Fig. 10. We arrive at 
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mu. 37202 



Determination of r ^ (r , E^,, 0p 0^ ) 



Fig. 14. 



-40- 



the relationship 
2 



^ = r ^ + E^ £l - sin^ 9^ cos^ (180 - 0^-2 ^E,^ [l - sin^ 9^ cos^ (180 - 0^)j 



- 1 : 



but from Fig. 14, cos to = 



E<p cos 9^ 



Erp j 1 - sin^ 9j cos 4 (180 - 0i)j 



7 TT 1/1 



substituting this value into the above relationship we find that 



r = /r(f + E^ fl - sin^ 9, cos^ (180 - 0-, )i - 2 r„ E_ cos 9, ( 

1 A/ u T ^ ^ J U1 1 



The possibility of a "skew 11 condition must likewise be considered 
in case of the Erj, vector. The length (^rp) must be determined and 
substituted into Eq.,16. The partial derivative of r^ with respect to 
.X(Erj,)must be taken and set equal to zero, 
thus : 

3 r l l v 2{% . 2 ~ 2 



= — 1 r 



S^(E t ) 2 [ 0 



t./?(E^)^ |1 -.sin~ 0^ cos^ (180 - 0 )] - 2 Tq ^(Erp) cos ( 



r 



[ 22 7 

1- sin 9^ cos (180 - 0^)j - 2 rg cos 9^ 

L 

Solving forj?^, we find that 



= 0 



Jl T = 



rg cos 9 



1 



["l - sin^ 9-^ co s ^ (180 - 0^)j 



(] 



for 0 < jlr£ < i(Erp) 



substituting this value back into Eq. 16 and calling this closest point 



of approach r^, we find that 
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1 Wo + ' 

A / 



or 



2 j:q cos^ [l - sin 2 9-] cos 2 (180- 0]_)J 2 r 2 cos 2 9 ^ 

E - sin^ 9i cos^ (180 - 0-, 

L 1 1 J 



|^1 - sin 9^ cos 2 ( 180-01 )j 



(18) 



3 = r o 1 " 
n 



cos 9i 



[l - sin 2 0! cos 2 (ISO - 0 X )J 



Summarizing the calculations for radial distances from the center 
wire, given E^, Ep, and r^ we must calculate four additional radial 
distances for each Ep, E^, vector pair. 



= t Q - Ep 



(9) 



= r. 



r^ + Ep ; JL - sin^ 9^ cos 2 (- ^ r o^P cos 



(ID 



1 - 



cos^ 9- 



' 1 - sin^ 9 2 cos 2 (- 02 /, 



(13) 



(only necessary to compute r A where — 5.^2^ 
and 0<-^p C Ep) 



X 1 



= / r^J + E 2 -[l - sin 2 9 1 cos 2 (180 - 0 X )] - 2 r Q E T cos 9 



2 /! 



(16) 



and; 




cos 9 



jE - sin^ 9^ cos' 



(ISO -0 x )j 



9 ^ = 180-92 — sin 



_]_ ; sin cos 2 T 



(18) 



. 2 . 2 . i/2 1 

ijl - sin <%,. sin a) j 

L L J 



where cosct= - 



(2E t tQ) 
2/3 E t Ep 



(14) 



